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We show that nano-scale variations of the order parameter in strongly-correlated systems can 
induce local spatial regions such as domain walls that exhibit electronic properties representative 
of a different, but nearby, part of the phase diagram. This is done by means of a Landau- 
Ginzburg analysis of a metallic ferromagnetic system near an antiferromagnetic phase boundary. 
The strong spin gradients at a wall between domains of different spin orientation drive the 
formation of a new type of domain wall, where the central core is an insulating antiferromagnet, 
and connects two metallic ferromagnetic domains. We calculate the charge transport properties 
of this wall, and find that its resistance is large enough to account for recent experimental results 
in colossal magnetoresistance materials. The technological implications of this finding for 
switchable magnetic media are discussed. 
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There is growing evidence that electronic spatial inhomogeneity on the nanoscale is an 
intrinsic property of strongly-correlated oxide materials. In one sense, this is not surprising, since 
the strong competing interactions produce many phases whose free energies are nearly 
degenerate, resulting in rich phase diagrams. Nevertheless, in particular cases, the origin of the 
phenomenon is often unclear, and its appearance is usually unpredictable and uncontrollable. 
The influence of this inhomogeneity on the macroscopic properties can be profound. 

A very important example is the family of high-Tc superconductors. The proximity of 
antiferromagnetic (AF) and conducting phases has led many authors to propose spatially 
inhomogeneous states, usually with a striped configuration [1]. For some materials and doping 
concentrations, there is unambiguous evidence for these inhomogeneities [2]. One may extend 
these ideas and propose that the cores of vortices, which are induced inhomogeneities, may be 
antiferromagnetic, i.e., may show signs of the proximity of an insulating phase in the doping- 
temperature phase diagram [3]. 

Another important class of oxides, the colossal magnetoresistive (CMR) materials, have very 
solid evidence for electronic inhomogeneity. Transmission electron microscopy has indicated [4] 
spatial charge separation, and scanning tunneling microscopy measurements [5] have provided 
strong evidence for a temperature range of local electronic phase separation even in nominally 
homogeneous ferromagnetic (FM) CMR materials. Measurements of electronic transport through 
magnetic domain walls in single-crystalline CMR materials [6] indicate a domain-wall resistance 
orders of magnitude greater than that attributable to scattering from the spin gradient in a Bloch 
wall, which is generally a very small effect [7,8]. We propose in this Letter that this phenomenon 
is the result of an AF core in the domain wall. Since the AF phase is insulating, the effect on 
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transport is dramatic. This has important technological implications. Domain-wall populations in 
films can be controlled with relatively small magnetic fields: they can be induced, as well as 
natural, inhomogeneities. Thus very large controllable magnetoresistances can be envisioned. 
This is a central requirement for magnetic memory elements. 

The bulk phase diagram in the temperature-doping plane of the CMR materials indicates a 
transition, usually first-order, from a metallic FM phase to an insulating AF phase. In the FM 
phase the order parameter (OP) has degenerate configurations, and domains of different ground 
states occur. In a domain wall, where the OP interpolates between the two configurations in a 
narrow region, there are strong gradients, and the bulk phase diagram no longer applies. 

If the system is close to a bulk phase boundary, it can be energetically favorable for a region 
of the wall to collapse into a different electronic phase rather than supporting a strongly rotating 
OP. Our calculations show that this can have drastic effects. It is possible to change the character 
of the domain wall from metallic to insulating, or even to dramatically change the spatial extent 
of the wall. As the system approaches the phase boundary, the width of the wall can diverge, or 
the walls could proliferate. This growth in spatial extent could be the origin of electronic phase 
separation observed in some of these materials. This would clearly affect macroscopic properties 
when the domain walls make up a significant fraction of the sample. Thus our theory of domain 
walls in films tends to confirm theories of bulk CMR which depend on the existence of 
mesoscopic patterns in which metallic and insulating regions interpenetrate [9]. 

Here we investigate the case of a single domain walls between FM domains in a material 
close to an AF phase boundary. However, the theoretical part of the paper is more general, and 
would apply to domain walls in other systems with closely competing OP's. 

II. Electronic phase of the domain wall 
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We model this system with a one-dimensional Landau-Ginzburg free energy 

F=-a^M' +b„M'-cMl + g„{M'f -a^+bA' + SaW +8^'^' (1) 
Here M is a vector OP for ferromagnetism, and A is a scalar OP representing the competing 
antiferromagnetism. Primes denote spatial derivatives. A is a scalar only for the sake of 
simplicity. The results for vector A would be very similar. We consider temperatures well 
below the FM ordering temperature, where all coefficients in Eq. (1) are positive. Furthermore, 
we are in a doping range where the bulk state is FM, but close to the I'^'-order boundary to the 
AF state: < / = F^*"' / ' < 1 > where [ F^"' ] are the bulk, uniform, free energy densities of 
the AF [F] states in the absence of the other. The anisotropy term -cM^ favors alignment of the 
ferromagnetic moment in the ±£ direction. This is the appropriate anisotropy for the PrSrMnOs 
films that we will be concerned with below. The coexistence term gA^M increases the free 
energy when the two order parameters are nonzero at the same spatial location. We will consider 
mainly the case where « • This is physically reasonable in the CMR systems since FM 
arises from the rather nonlocal double exchange mechanism. 

We impose a magnetic domain wall by enforcing the boundary conditions M — > ±M^"^z at 
X = +00 . The wall itself itself is only thermodjoiamically stable in the presence of long-range 
(dipole or other) forces [9]. Here M*"' =((a„ + c)/2Z7^)"' is the equilibrium FM magnetization 
of the uniform, single-domain system. We explore broad regions of parameter space using 
simulated annealing of the free energy in Eq. (1), and present analytic solutions for particular 
regions of phase space. 

With no coupling of the two OP's (g = ), Eq. (1) describes a ferromagnetic Bloch wall. Then 
gm defines the spatial length scale, and bm becomes a scale factor for the overall free energy. The 
qualitative characteristics of the wall are controlled only by am- For am large the magnitude of the 
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magnetization vector is stiff, and Eq. (1) can be analytically minimized to obtain the known 
result for a Bloch domain wall where M rotates at constant magnitude. This gives 
^ = arctan(exp(x/^), where is the angle of the local magnetization with respect to the 
quantization axis, and ^= (Zg^/cY'^ . For smaller am, the magnitude M is reduced at the wall 

center, and also M completes its rotation over a slightly shorter distance. 

Nonzero g results in competition between the two order parameters. In this case we find an 
additional solution for large enough g, shown in Fig. 1 (for = and g=S), where an 
antiferromagnetic region (dashed line) exists in the core of the domain wall. In this case the 
ferromagnetic magnetization M (solid lines) goes continuously through zero. The general 
characteristics of this type of solution are robust over a wide parameter range: M does not rotate 
continuously but rather is directed along ±£, changing orientation discontinuously as the 
magnitude passes through zero. Solutions for two different parameter sets are shown: the 
'narrow' wall corresponds to parameters such that/ = F^*"' /Fp"^ = 0.64, (relatively far from the 
bulk phase boundary), while the 'wide' wall has/=i^"^ /Fp"^ = 0.98 (close to the boundary). The 
absence of an antiferromagnetic gradient term (g^ = 0) leads to a discontinuous slope as A goes 
to zero. This sharp change allows a straightforward definition of the antiferromagnetic core 
width as the region over which A is non-zero. Figure 2 shows that this width diverges 
logarithmically as / — > 1 . 

These solutions have = and = m . In this case, the free energy can be expressed in 
normalized variables to clarify the physics. Equation (1) becomes 

F = [-M' + M*/2 + (M'f /2] + f'\-A^ + A*/2] + ')A^M\ (2) 



5 



where M=MJM'"\ A = A/ A'"> , f = / p;;' , F=FI{Fl^'Fr) , ^ = g/(4Z7^^)"% and the 
normalized spatial coordinate x = x/^^ with = (2g^/(a^ +c))"^. The shows that the free 
energy, and hence the physical behavior, is determined by only two parameters: the ratio of the 
free energies for the isolated uniform solutions and the coexistence energy. 

This form permits an analytic calculation of the domain wall structure. Making Eq. (2) 
stationary with respect to variations in A and M leads to 

A' = l-(r//"2)M2 

(3) 

M'72 = - M(1-/"V') 
Equation (3) shows that A decreases as increases from zero and, for If large 
enough, A = beyond a critical ferromagnetic magnetization value m = f'^lj- This divides 
the solutions to Eq. (3) into one valid in the 'exterior' (A = 0) region and one valid in the 
"interior" (A^^^O). In the exterior region, Eq. (3b) can be integrated (using the boundary 
condition M' — > as M — > 1) to find M^^ = tanh(x - . x^ is determined by matching to the 
interior solution. The interior solution can be integrated analytically ( after substituting for A and 
matching the interior and exterior solution derivatives at m = f'^ If ) to obtain 

{M'f ={\- M'f - f{\-{rl f")M'f , 
which can then be integrated numerically to find x{M) as 



dM 



X 



= \ , (4) 



The width of the domain wall antiferromagnetic core is obtained by integrating up to 
M= if^'^ I yi'^ , where A goes to zero, and multiplying by two. For large y the wall width 
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diverges as / — > 1 , which in Eq. (4) arises from the divergence of the integrand at M=0, 
/ — > 1 . The leading-order behavior of Eq. (4), obtained by approximating the integrand for small 

M, is -(l/«)ln(l-/) + (2/«)ln«, with a= -^lif^'^/ -I) . This reproduces the logarithmic 
dependence on 1 - / , and the ~ ^"^'^ dependence of the slope, shown in the inset to Fig. 2. 

From a simulated annealing analysis, we find that over a broad parameter range this 
antiferromagnetic core solution is very close to the one that minimizes the general free energy 
(Eq. (1)). This is illustrated by the 'domain wall phase diagram' of Fig. 3, where the minimum 
free energy solutions obtained by simulated annealing of Eq. (1) are characterized for different 
values of the parameters / and y . The region labeled 'AF-core' is a region of parameter space 
where the antiferromagnetic OP A is nonzero in the domain wall. Solutions in the 'F-core' region 
have A = throughout the wall (Bloch wall solution). The solid lines are first-order phase 
boundaries, based on the numerically observed metastability (hysteretic behavior) of the two 
distinct solutions on either side of the boundary. The dashed line is a second-order phase 
boundary, again based on simulated annealing results. Also shown in Fig. 3 are phase boundaries 
for various values of the antiferromagnetic gradient coefficient g^, indicating that the generic 
behavior is preserved away from the special case = 0. 

The AF cores in CMR materials, if they exist, would be insulating, as in the bulk. This would 
lead to high domain wall resistances, as is in fact observed. This resistance can be calculated in 
simple model of the band structure. Our model will be a classical spin S to represent the 3 
localized electrons on each Mn site, and a single tight-binding band to represent the current- 
carrying delocalized eig electrons. The Hamiltonian for the electrons is then 

Here / is a site label and is the spin of the conduction electron. 
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J is the Hund's rule coupling. The spins 5, of the electrons are held fixed. We adopt a 

picture of the domain wall where S^ = M(x)z + A(x)x, and the position dependence of M and A 
is given in Fig. 1. A Neel structure is assumed for the cubic lattice of Mn atoms. 

We compute the spin-dependent (multichannel) transmission matrix numerically. From this, 
the resistance of the domain wall follows by using a generalized Landauer formula. The 
resistance as a function of wall width is shown in Fig. 4. The resistance is roughly exponential. 
This can be understood within a WKB tunneling picture. If the Fermi level Ef of the metal lies 
in the gap of the insulator, then the wavefunction ^is damped in the wall: y/ ~ exp(x/ x^) in the 

wall, where x^ = at/\EQ- Ep\. Here -£^^1 is the smallest distance from the metallic Fermi 

level to a gap edge in the insulator. 

The calculation presented here gives the conductance per unit area of a single domain 
wall. In experimental measurements [6], the excess electrical resistance of a sample with 
multiple magnetic domain walls is compared to that of the sample in the single-domain state. The 
ratio of these (F) is related through simple geometrical effects to the ferromagnetic domain 
width d and the single-domain conductivity <T^ as F = <T^ /g^d. Thin films of Pro.eTSro.asMnOa 
experimentally show large domain wall resistivity [6], are strained by the substrate, with a Curie 
temperature lower than unstrained material. Strain has moved the sample closer to the 
antiferromagnetic phase boundary. In addition, the strain generates an out-of-plane magnetic 
anisotropy and corresponding dense domain structure. In this case, the ferromagnetic domain 

o 

size is on the order of the film thickness, typically 200 A in these thin films. 

Experimentally the ratio F is found to vary with temperature, but can be very large [6]. We 
can get an approximate upper bound for the wall conductance per unit area by noting that 
g^ = (e^/h)n^'^{T), where n = 5.56xlO^Vcm^ is the electron density and {T)<1 is the 
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appropriately averaged transmission coefficient. If we take a domain conductivity 
£7^ =10'*/(^2-cm) typical of manganite systems and d = 200 A, we find r>0.4. Thus the wall 
resistance is comparable to the domain resistance at small thicknesses, and Fig. 4 shows that the 
wall resistance can easily dominate the entire resistance of the film. This is due to the presence of 
a switchable metal-insulator transition in the boundary between two conductors. 

The presence of local regions of insulator on the metallic side of the transition also suggests 
some unusual routes for the first-order transition not only in thin films but also in bulk samples. 
Let us focus first on the case of tetragonal symmetry, for then the topological character of the 
problem simplifies: there is only one type of wall, which separates only two possible equivalent 
domains. A standard scenario as the transition is approached from the ferromagnetic side would 
be nucleation of nearly spherical regions of insulator in a homogeneous metallic phase. 
However, the ferromagnet will always have domains. (Even in equilibrium, the dipole interaction 
enforces this.) As the system approaches the transition, the new phase could nucleate in the 
walls; the insulating regions would be topologically planar rather than spherical, and would 
traverse the entire sample. Two routes through the transition can now be envisioned. As their 
local free energy is lowered, (1) the walls could proliferate and eventually join; (2) the wall 
widths could increase continuously (as suggested by Fig. 1) resulting in coalescence. Scenario 
(2) would be the result if the only forces in the problem are those represented by Eq. (1) plus 
dipole forces. However, higher-order gradient and anisotropy terms can set a fixed length scale 
for the wall width [10], which would promote scenario (1). 

If the S5rametry is higher, say cubic, then several different t5^es of walls are possible and 
complex networks may drive the transition [9]. 
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FIGURE CAPTIONS 



Figure 1: Typical solutions for the spatial dependence of the ferromagnetic (M) and 
antiferromagnetic (A) order parameters at a wall between oppositely oriented ferromagnetic 
domains. The parameter / indicates the proximity of the material to the antiferro-ferro phase 
boundary in the bulk system. 

Figure 2: Width of the antiferromagnetic domain- wall core as a function of proximity to the 
antiferro-ferro phase boundary for 7 = 2. The inset indicates the logarithmic divergence for 
different values of the coexistence coefficient y . The dashed portion of the curve indicates the 
regions where the ferromagnetic domain wall solution is of lower free energy. 

Figure 3: Phase diagram of the system, showing regions in which the domain wall has an 
antiferromagnetic, insulating core. The three phase boundaries shown are for different values of 
the antiferromagnetic gradient coefficient. 

Figure 4: Interface resistance of the domain wall as a function of wall width. The width is 
defined as the length over which the antiferromagnetic moment A is nonzero, measured in unit 
cells. The resistance grows roughly exponentially, with minor interference effects superposed. 
The calculations are for a tight-binding model with an electron density of n = 5.56xlO^Vcm^, 
/ = 0.88 eV, and t = 0.27 eV as estimated from band-structure calculations [11]. 
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Figure 1: 
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Figure 3: 
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Figure 4: 
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